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Abstract 

We analyse various flavour changing processes like t —)> hu, he, h —>■ re, rg 
as well as hadronic decays h ^ bs, bd, in the framework of a class of two Higgs 
doublet models where there are flavour changing neutral scalar currents at tree 
level. These models have the remarkable feature of having these flavour-violating 
couplings entirely determined by the CKM and PMNS matrices as well as tan /3. 
The flavour structure of these scalar currents results from a symmetry of the La- 
grangian and therefore it is natural and stable under the renormalization group. 
We show that in some of the models the rates of the above flavour changing 
processes can reach the discovery level at the LHC at 13 TeV even taking into 
account the stringent bounds on low energy processes, in particular g ey. 
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1 Introduction 


The second run of the LHC, at a center of mass energy ^/s = 13 TeV, will provide an 
important probe of flavour changing couplings of the recently discovered scalar boson 
h gig. These couplings can contribute to rare top decays like t ^ hq (q = u, c) 
and may also lead to leptonic flavour changing decays such as h —)■ (£ = p, e), 

as well as hadronic flavour-changing decays like h ^ bs and h —)■ bd. These decays 
are strongly suppressed in the Standard Model (SM) since these couplings vanish at 
tree level. However, Higgs Flavour Violating Neutral Couplings (HFVNC) can arise 
in many extensions of the SM, including in Two Higgs Doublet Models (2HDM) gg. 
Any extension of the SM with HFVNC has to comply with the strict experimental 
limits on processes mediated by flavour changing neutral currents as well as with the 
limits on CP violating transitions leading, for example, to electric dipole moments of 
quarks and leptons g. 

In this paper, we investigate the allowed strength of HFVNC in the framework of a 
class of 2HDM, denoted BGL models, hrst proposed for the quark sector g, generalised 
in g and then extended to the leptonic sector g. These models have the remarkable 
feature of having HFVNC, but with their flavour structure entirely determined by 
the Cabibbo-Kobayashi-Maskawa (CKM) j9,10 and the Pontecorvo-Maki-Nakagawa- 
Sakata (PMNS) 01 ig matrices, denoted V and U respectively. HFVNC have been 
widely addressed in the literature 13-41 . The distinctive feature of BGL models 


is the fact that many of the most dangerous HFVNC are naturally suppressed by 
combinations of small mixing matrix elements and/or light fermion masses. This is 
achieved through the introduction of a symmetry in the Lagrangian and therefore 
the suppression is entirely natural. Another salient feature of BGL models is that 
depending on the specihc model of this class, HFVNC exist either in the up or in the 
down sector, but not in both sectors simultaneously. Analogous considerations apply 
to the leptonic sector. We will pay special attention to the evaluation of the maximum 
rate at which these processes can occur in this class of models, without violating the 
stringent bounds on processes like /i —>■ ey. In the general 2HDM, in the so-called 
Higgs basis 42-45 , there are three neutral scalars, which we denote H°, R° and A. 
The couplings of H*^ to fermions are flavour diagonal in the fermion mass eigenstate 
basis. On the other hand, in the general 2HDM the fields R° and A have HFVNC with 
arbitrary flavour structure. The remarkable feature of BGL models is the fact that the 
flavour structure of HFVNC only depends on V and U. Furthermore, the neutral scalar 
mass eigenstates are linear combinations of H° and R° together with the pseudoscalar 
neutral held A. In these models the imposed symmetry restricts the scalar potential 
in such a way that it cannot violate CP either explicitly or spontaneously, therefore 
once we perform a rotation through the angle fd that takes the helds from the basis 
chosen by the symmetry to the Higgs basis, the charged helds and the pseudoscalar 
held A are already physical helds. As a result, the two other neutral physical helds 
are related to H° and R° through a single angle rotation. In a previous work 


46 


we have performed a detailed analysis of the allowed mass ranges for the new scalars 
under the assumption that the discovered Higgs h can be identihed with H°. We have 
then shown that in some of the BGL models these masses can be in the range of a few 
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hundred GeV and thus within reach of the LHC 13 TeV run. In this paper, we work 
in the general case where h is a mixture of H° and R°, controlled by an angle denoted 
(3 — a. The strength of the HFVNC of h depend crucially on tan/3 = with n* 

the vacuum expectation values of the scalar doublets, and cos(/3 — a). BGL models 
have many features in common with several of the implementations of the minimal 
flavour violation hypothesis (MFV) (W ,47 49 . However, BGL models have the unique 


feature of coming from a symmetry, and as a result they have a reduced number of free 
parameters. This allows for dehnite predictions once constraints on these parameters 
are imposed, taking into account the present experimental bounds. 

The challenge is to answer the following question: in some of the BGL models, can 
one have regions, in the tan/3 versus a — 13 plane, where the HFVNC of h are such 
that the rare processes t —)■ hq, h ^ fir can occur at a rate consistent with discovery 
at LHC-13 TeV? Of course, these regions have to be consistent with the stringent 
constraints on all Standard Model (SM) processes associated to the Higgs production 
and its subsequent decays in the channels ZZ, WW, 77 , bb and rr. Furthermore the 
constraints derived from low energy phenomenology have to be considered: both those 
obtained in 46 and those new due to the presence of H°-R° mixing. Processes such as 


h ^ bs and h ^ bd are probably not within reach of the LHC but become important 
for the physics of the future Linear Collider. In this paper, we also address the question 
of what BGL models could lead to the observation of such decays by the future Linear 
Collider. 

The paper is organised as follows. In the next section we briefly review the main 
features of BGL models in order to set the notation. In Section we analyse top 
flavour changing decays of the type t ^ hq {q = u, c) in the framework of BGL models 
with HFVNC in the up sector. In Section]^ we perform an analysis of flavour changing 
Higgs decays in BGL models. In particular, we consider neutrino models with HFVNC 
in the charged lepton sector giving rise to h —)■ f'r decays. Up type models are also 
considered, giving rise to h ^ bs flavour violating decays. In Section we investigate 
the discovery regions and the existing correlations for these decays in the framework 
of BGL models. In section we present our conclusions. The paper contains two 
appendices, in Appendix we explain how the relevant Higgs experimental data has 
been incorporated into the analysis, and in appendix B.3 we give details relative to 
the low energy flavour constraint /r —>■ ey. 
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2 Main Features of BGL Models 


Our work is done in the context of Two Higgs Doublet Models. The quark Yukawa 
interactions are denoted by: 

Cy = -Wl [ri + Ts $2] dl-WL [Ai <i>i + A21>2] < 

- Zf [Hi $1 + n 2 $ 2 ] - Zf [Si $1 + S 2 $ 2 ] 4 + h-c-, 

where Tj, A* Hj and Sj are matrices in flavour space. The requirement that Tj, Aj lead 
to tree level Flavour Changing Neutral Couplings (FCNC) with strength completely 
controlled by the Cabibbo - Kobayashi - Maskawa matrix V, was achieved by Branco, 
Grimus and Lavoura (BGL) by means of a symmetry imposed on the quark and 
scalar sector of the Lagrangian of the form: 

Qlj ^ exp (ir) Qlj , ^ exp (i 2 r) , $2 ^ exp (zr) $2 , ( 2 ) 

where r 7 ^ 0 , vr, with all other quark helds transforming trivially under the symmetry. 
The index j can be hxed as either 1, 2 or 3. Alternatively the symmetry may be chosen 
as: 

Qlj ^ exp (ir) Qlj , ^ exp (i 2 r) d%j , $2 ^ exp (-ir) $2 . (3) 

The symmetry given by Eq. ( p|) leads to Higgs FCNC in the down sector, whereas 
the symmetry specihed by EqZ^ leads to Higgs FCNC in the up sector. These two 
alternative choices of symmetry combined with the three possible ways of hxing the 
index j give rise to six different realizations of 2HDM with the flavour structure, in 
the quark sector, controlled by the V matrix. We call up-type BGL models those 
with HFVNC in the down sector, coming from the symmetry given by Eq. ([^ and we 
identify each one of the three implementations by u-type, c-type or t-type depending on 
the value of the index j, respectively 1, 2 or 3. Likewise for the down-type models. In 
the leptonic sector with Dirac neutrinos there is perfect analogy with the quark sector 
and the corresponding symmetry applied to the leptonic helds leads to six different 
realizations with the strength of Higgs mediated Favour changing neutral currents 
controlled by the Pontecorvo - Maki - Nakagawa - Sakata matrix, U. As a result there 
are thirty six different implementations of BGL models. As was shown in reference [^, 
in the case of Majorana neutrino there are only 18 models corresponding to the neutrino 
types and therefore with HFVNC in the charged lepton sector. 

The discrete symmetry of Eqs. ([^ or (|^ constrains the Higgs potential to be of the 
form: 

V = + /i 2 <h ^$2 - mu ($ld >2 + $ 2 *^* 1 ) + 2 A 3 j (^$^<1)2 

+ 2A4 ($I<I>2) + Ai + A2 ($^$2)', 

the term in mi 2 is a soft symmetry breaking term. Its introduction prevents the 
appearance of an would-be Goldstone boson due to an accidental continuous global 
symmetry of the potential, which arises when the BGL symmetry is exact. Such a 
potential cannot violate CP either explicitly or spontaneously. As a result the scalar 
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and pseudoscalar neutral Higgs fields do not mix among themselves and we are left 
with only two important rotation angles, (3 and a. The angle f3 is such that the 
orthogonal matrix parametrised by this angle leads to the Higgs basis, singling out the 
three neutral fields: H°, with couplings to the quarks proportional to mass matrices, R° 
which is a scalar neutral Higgs and A which is a pseudoscalar neutral Higgs; as well as 
the physical charged Higgs fields and the pseudo-Goldstone bosons. In BGL models 
the fields A and are already physical Higgs fields, while H*^ and may still mix. 
In the limit in which H° does not mix with R°, H*^ is identified with the Higgs field h 
recently discovered by ATLAS Q and GMS . In this limit this field does not mediate 
tree level flavour changes and a is defined in such a way that the mixing angle between 
these fields, {(3 — a), acquires the value 7t/2. In fact expanding the neutral scalar fields 
around their vacuum expectation values 50 in the form 0° = -^(vj + Pj +iT]j) we can 
write: 


H°\ _ / cos (3 sin (3 
R° j ~ I—sin/3 cos/3 


cos a sm a \ l pi 
-sina cosay \p 2 


( 5 ) 


The angle /3 is determined by tan/3 = V 2 /V 1 ] in the following we use the shorthand 
notation tan/3 = cos(/3 — a) = C/ja and sin(/3 — a) = spa- 

The general form for the Yukawa couplings of 2HDM written in terms of quark 
mass eigenstates and the scalar fields in the Higgs basis is given by: 


— — 


y/2H^ 


u 


- NIv^l 


d + h.c. 


_H0 

V 

_ RO 

V 

A 

+ 

V 


\u Du u + d Ddd] 

u {Nu7r + u + d (A^rf7H + ^]7 l) d 

u (A'„7r - Nl-iL) u-d - Nl^iL) d 


( 6 ) 


where 7 ^ and 'Jr are the left-handed and right-handed chirality projectors, respectively, 
and Dd and Du are the diagonal mass matrices for down and up quarks respectively. 
This equation defines the matrices and which give the strength and the flavour 
structure of FGNG and are also involved in the couplings of the charged Higgs fields. In 
general 2HDM, still in the Higgs basis, the flavour structure of the quark sector is fully 
determined by the quark masses, the V matrix and the two matrices and N^. It is 
worth emphasising the high predictive power of the general 2HDM, as can be seen from 
Eq. 0. Let us assume that a pair of charged Higgs and the three neutral scalars (in 
general the physical neutral scalars are combinations of H°, R'^ and A) are discovered. 
The couplings of to quarks can be readily measured from their decays. Since V 
in Eq. 0 stands for the GKM matrix which is known, from decays one can derive 
and N^, which enables one to predict in the framework of the general 2HDM the 
flavour structure of the neutral scalar couplings. This would be essential to prove that 
the discovered neutral and charged scalar particles were part of a two Higgs doublet 
structure. In general 2HDM, the matrices and N^, are entirely arbitrary. On the 
contrary, BGL models have the remarkable feature of having the flavour structure of 
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and entirely determined by fermion masses, V and the angle j3 with no other free 
parameters. This results from the symmetry introduced in the Lagrangian, in order to 
achieve the BGL flavour structure in a natural way. As previously emphasized, each 
one of the six implementations of BGL in the quark sector only has FGNG in one of 
the quark sectors, either up or down. In BGL up-type models the matrices and 
have the following simple form: 

= [*/><*••. - + tf'KrVi.] (D,)„ , (7) 

where no sum in j implied. The upper index [uj) indicates that we are considering 
a symmetry of the form given by Eq. ([^, i.e., an up-type model with index j, thus 
leading to FGNG in the down-sector. Notice that all FGNG are proportional to the 
factor {tp + multiplying products of entries involving one single row of V. The 
corresponding matrix is given by: 


is a diagonal matrix, the tp dependence is not the same for each diagonal entry. 
It is proportional to for the (jj) element and to for all other elements. The 
index j Exes the row of the V matrix which suppresses the flavour changing neutral 
currents. Since, for each up-type BGL model a single row of V participates in these 
couplings, one may choose a phase convention where all elements of and are 
real. For down-type models, which correspond to the symmetry given by Eq. ([^, the 
matrices and exchange role, and now we have: 

(N'-T, = [HK. - (H + t-f'KiV,)] (Dj„, (9) 

= [«« - (‘s + (£>.)„«„. (10) 

In down-type models the flavour changing neutral currents are suppressed by the 
columns of the V matrix. 

In this paper we allow for the possibility of h being a linear combination of H° and 
R° which is parametrised by the angle (/3 — a): 

h = + C^aR° , ff = C^aH° - S^aR° ■ (H) 

This mixing is constrained by data from the LHG Higgs observables (see appendix [A]) . 
The quark Yukawa couplings of the Higgs field h can be denoted as: 


^hqq ^ij ^Li ^Rj h ^ij ^Li '^Rj T h.C. , 


( 12 ) 


and similarly for the leptonic sector with the coefficients denoted by Yfj and From 
Eqs. (|^ and ([II|), we can read: 


Y.^ = - 

V L 


T (A(. 


dJV 


Y^j — \_S/3a {Du)ij + C/3q. {Nu)ij~\ 


(13) 


More specihcally, for i ^ j, we get the following flavour violating Yukawa couplings, 
for the different types of BGL models: 


5 




(i) up-type Uk model, with k fixed as 1 {u) or 2 (c) or 3 (t), where HFVNC arise in 
the down quark sector: 

^ Cf)a{tp + , z 7 ^ j, no sum in k , (14) 

(ii) down-type dk model, with k fixed as 1 (d) or 2 (s) or 3 ( 6 ), where HFVNC arise 

in the up quark sector: 

TTL 

^ij ( 4 ) = -VikV*k + t-p^) , z 7 ^ j, no sum in k , (15) 

(iii) leptonic sector, neutrino-type, model, with k fixed as 1 {vi) or 2 (z/ 2 ) or 3 ( 1 / 3 ), 

where HFVNC arise in the charged lepton sector: 

TTli 

YlMk) = -UikUjk , z 7 ^ j, no sum in k . (16) 


In the case of Dirac neutrinos one can write similar expressions for charged lepton type 
models and in this case the FCNC appear in the neutrino sector and are suppressed 
by the extremely small neutrino masses. In the case of models of the charged lepton 
type, only diagonal couplings are relevant. These couplings, as all other diagonal ones, 
can be extracted from equations ([^ - ( 10 ), and if necessary making the corresponding 
changes of quarks by leptons. 
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3 Flavour changing decays of top quarks 


In this section, we analyse flavour changing decays of top quarks t — hq. They can arise 
in down-type BGL models, where there are Higgs flavour violating neutral currents in 
the up sector. According to Eqs. (13) and (15), the couplings of the Higgs particle h 
with a top t and a light up-type quark u or c, in a model of type dp, can be written as 



+ q = li,c. 


(17) 


One can then evaluate the corresponding t ^ hq decay rate. As previously mentioned, 
there are three types of models of this class, dp, depending on the column of the V 
matrix which suppresses the flavour changing currents. The result is 


^idp)it -)■ hq) 


mf 

3271v‘^ 





2„2 


cgJh + ts ) 


(18) 


Note that, apart from the global factor +1^) , every other factor in Eq. (18) is 

fixed once we choose the specific down-type model dp and the decay channel t ^ he 

or a given model, t ^ hq processes constrain the factor 
we enumerate the decay channels and the models according 




hu. Therefore 

-1x2 


or t —)■ 

.2 

''pOL 

to the V factors involveoT 


In Table 


Model 

t ^ hu 

t ^ he 

d 

in,!). 

2 A®) = 7.51 ■ 10-5 

\V^V„ 

2 A^) = 4.01 ■ 10-6 

s 

\V^,Vu 

2 A®) = 8.20 ■ 10-® 

\V,,V,. 

2 A^) = 1.53 ■ 10-3 

b 


2 A^) = 1.40 ■ 10“^ 


2 (A. A^) = 1.68-10-3 


Table 1: V factors entering Eq. ( |I^ , A ~ 0.22 is the Cabibbo angle or Wolfen- 
stein main expansion parameter |52|. 


It is clear that the most interesting models for t ^ he are the s and b models, where 
the suppression is only at the level, compared to the d model which has a strong 
suppression for the same decay at the A® level. The d model has the curiosity that the 
suppression is higher for t ^ he than for t —)■ hu, unlike in the other two models. The 
branching ratio for f —)■ hg in the dp type model is 


Br(d,)(t ^ hq) = 


t hq) 


Wb) 


— f{xh,yw) 


+^/ 3 ^)^ ( 19 ) 




tb I 


where 


1 -1 

f{xh, yw) = 0(1- ^hf (1 - ^yw + 2|/w) > with Xh = yw = ■ (20) 


Using the top quark pole mass rrit = 173.3 GeV 51 , rrih = 125.0 GeV and = 80.385 
GeV, one obtains f{xh,yw) = 0.1306 . Gonsidering then the upper bounds 0.79% from 
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the ATLAS 53 and 0.56% from the CMS 54,55 collaborations, we obtain, for 6 and 
s-type models, the following constraint 

|c/3Q(i/3 + < 4.9 . (21) 

Notice that for this valne, pertnrbative nnitarity constraints have to be considered (see 
appendix |B.1[) . 


4 Flavour changing Higgs decays 

4.1 The decays h ^ ir {i = jji, e) 

The most interesting BGL models with HFVNC in the leptonic sector are the u models, 
where there are FCNC in the charged lepton sector. As seen in the previous section for 
the quark sector, there are three neutrino-type BGL models, depending on the column 
of the U matrix which enters the FCNC in the leptonic sector. Using a notation 
analogous to the one of the quark sector and considering Eq. (16) for the h coupling to 
p and r, we have 




( 22 ) 


and the decay rate: 


/xr) + F(^^)(h ^ pr) = Tsuih^ rr), (23) 

2 

with FsM(h —>• rr) = Notice, again, the appearance of the same factor 




{tjS + t 


/3 


type models. 


Table lists the PMNS mixing matrix factors for the different v - 


Model 

h — )■ e/i 

h ^ er 

h ^ HT 


IC'elC'al 

^(~i) = 0.105 


^(~ i) = 0.118 

lUM 

^(~ ^) = 0.028 



2(^ 1) = 0.089 

|C'.2£'.2 

^(~ i) = 0.126 

|C7„2V2 

^(~ i) = 0.115 


WaU 

„s|" = 0.0128 

= 00097 

|f7„3Vs 

^(~ i) = 0.234 


Table 2: U factors entering Eq. (23) for the different u - type models; estimates, e.g. 
1/9, 1/36, correspond to a tri-bimaximal U (except, of course, for lUgsl); analogous 
information for h —)■ e/r and h ^ er decays is provided. 


The hrst direct search for lepton-flavour-violating decays of the observed Higgs 
boson performed by the CMS collaboration 56 , led to the observation of a slight 


excess of signal events with a signihcance of 2.4 standard deviations. The best ht value 

is: 

Br(h —)■ pr-F r/2) = (O.Sdlgl®) % , (24) 


which sets a constraint on the branching fraction Br(h —)■ fir + Tfi) < 1.51% at the 95% 
conhdence level. The ATLAS collaboration has presented a result based on hadronic 
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r decays giving Br(/i —)■ /xr + r/x) = (0.77 ± 0.62)%. Assuming the h width to be 




[SM]^ 4.03 MeV), one can use the SM branching ratio BrsM(h —)■ rr) = 0.0637 
3), and obtain the estimate 

1, (25) 


in Eq. (23), and obtain the estimate 

\Cjiaitp + ) I 


necessary to produce Br(h —/ir + rp) of order 10 


-2 


4.2 The flavour changing decays h ^ bq {q = s, d) 


We now address up-type BGL models, where there are scalar mediated FCNC in the 
down sector and the most promising experimental signatures correspond to h ^ bq 
decays, with q = s,d. The relevant flavour changing h couplings to the down quarks 


in Eq. (12) are, according to Eq. (14): 


rrib 


Y,b (Uk) = -Cfaitg + tg ) 14', 14, — , qgtb, no sum in k 


( 26 ) 


Once again, it should be emphasised that once the up-type model Uk is chosen, the 
strength of the flavour changing couplings only depends on the combination Ci 3 a{ti 3 +t'^^) 
together with the down quark masses and V factors which are already known. The 
decay rate of h to pairs of quarks qiqj {i % j) is 


rK)(h ^ qiqj + qiqj) 


Stt 


1 

2 


yo\" + 


1 

2 



(27) 


We thus have 

-^bq + bq) = + t^^) Tsuih ^ bb). (28) 

Assuming that Th — we can make the following estimate 

BrK)(h ^bq + bq) = f \Vkq?\V^,^,\^ BrsM(h ^ bb) , (29) 

where BrsM(h —)■ bb) = 0.578. The relevant CKM-related factors for h ^ bs and 
h —)■ in all three Uk BGL models are given in Table 


Model 

h ^ bd 

h ^ bs 

u 

IK.,V4, 

|2(^ A®) = 1.33- lO-'^ 

IK., Vs 

2(^ A*) = 7.14-10-^ 

c 

IVsK* 

2 (^ A®) = 8.52 ■ 10-5 

IK.Ks 

2 A^) = 1.59 ■ 10-3 

t 

KgKb 

2 (^ A®) = 7.90 ■ 10-5 

IK.Ks 

2 (^ A^) = 1.61 ■ 10-3 


Table 3: V factors entering Eq. (28), A ~ 0.22. 


We thus have, to a good approximation: 

• in models c and t, 

Br(h -^bs + bs) ~ cJo(t/3 + t^^)^ ~ 10"^ + i^^)^ , (30) 
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• in model m, 


Br(/i6s + 6s) ~ + ~ 10 , (31) 

• in all M, c and t models, 

BT{h -^bd + bd) ~ A® ~ 10"® c^a(^/3 + • (32) 


We stress that, a priori, in models where there is no h —)■ pr constraint, one can reach 
values for Br(h —)■ 6s + s6) not far from 10“^. This can happen in charged lepton 
models of the charm and top types with cpaitp + ^^^) ranging from 5 to 10. Again, see 
appendix B.l for perturbative unitarity constraints on these values of 


5 Correlations among Observables 


One of the most interesting aspects of flavour violation in BGL models is the fact that, 
in this framework, it is possible to establish clear correlations between various flavour 
violating processes. As previously emphasized, one of the key features of the BGL 
models analysed in this paper is the presence of flavour changing neutral currents at tree 
level, but naturally suppressed by entries of the CKM and/or PMNS mixing matrices. 
Apart from these mixing matrices, in these models FGNG only depend on the values 
of tan /d, cos(/d — a) and fermion masses. However, the level of suppression depends on 
the specihc BGL model, which in turn implies that the correlations differ from model 
to model. In this paper we analyse the tree level flavour violating decays involving 
the Higgs boson already discovered at the LHG, which were listed in the previous 
section. It should be pointed out that the analysis has to take into consideration 
the flavour conserving Higgs constraints already obtained from Run 1 of the LHG. In 
particular one has to comply with the measured signal strengths for the hve relevant 
decay channels measured, to wit: h —)■ W^W~, h —)■ ZZ, h ^ bb, h ^ rf, and 
h — ?■ 77. They involve the flavour conserving couplings of the Higgs and put constraints 
in the vs. a — (3 available space. A detailed explanation of observables and the 
constraints is given in appendix An extended analysis of the phenomenology of the 
models under consideration, of the type presented in 46 , is beyond the scope of this 
paper. However we also take into consideration the current bounds from the low energy 
processes Bd — Bd mixing, Bg — Bg mixing, mixing and mixing. The 

experimental limits on these processes can be, in principle, easily translated into limits 
on the combination which is relevant for our models, as will be shown in 

what follows. As pointed out before, there are several different types of BGL models 
and depending on the model under consideration we have FCNC in the down sector 
or in the up sector, but never in both sectors at the same time. We can thus analyse 
different kinds of correlations among the different HFVNG observables considered in 
the previous sections: 


within the same quark sector: 
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— t ^ he vs. t ^ hu, in down-type models, where there are tree level FCNC 
in the up quark sector, 

— h ^ bs vs. h —)■ bd, in up-type models, where there are tree level FCNC in 
the down quark sector, 

• within the quark and the lepton sector in neutrino-type models (where there are 
tree level FCNC in the charged lepton sector), 

— t —)■ hg vs. h —)■ /ir, in down-neutrino-type models, 

— h ^ bqvs. h —)■ pr, in up-neutrino-type models. 


56 from the CMS collaboration 


57 - pointing towards the possible 


Taking into account the different features of specific BGL models, we will present 
various correlations in the framework of BGL models where the above processes may 
occur at a level consistent with discovery at the LHG and/or at the future Linear 
Gollider. Notice that, at present, there is evidence 
- not challenged by the recent ATLAS result 
observation of the decay h —)• /ir, which would constitute a clear signal of physics 
beyond the SM. If such an evidence persists it will favour neutrino-type BGL models. 
Top quark decays such as t ^ he and t ^ hu will be further explored in Run 2 of the 
LHG. Goncerning h ^ bs and h —)■ bd, these decays are probably inaccessible at the 
LHG due to the overwhelming backgrounds. However, they constitute very promising 
channels for the future Linear Gollider 
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5.1 


Correlations in Down - Charged lepton type models: t 
versus t —hu 


he 


We have chosen down - charged lepton type models in order to have FGNG in the up 
sector while avoiding flavour violation in the charged lepton sector. In the down-type 
BGL model d^, with HFVNG in the up sector mediating rare top decays, following 


Eq. (19), we get the simple correlation formulas 


Br(rf^)(f ^ he) _ 117^1 


Br(d^)(f^hM) | 1 /. 12 


(33) 


u'y \ 


Values of the ratio in Eq. (33) for different models are shown in table 


Model d.y 

d 

s 

b 

hvP 

/ , \2'| _ Q nc,o/| 

IRsP / , 1 \ _ 1 Q acA 

iRbh / , 1 _ 11 f) KQ 


\y 12 t ^ ) U.UOO^ 

1 ^ ud 1 


lyjA - a2(pW)^ 


Table 4: Ratio Br(rf^)(t —)■ hc)/Br^)(t — hu) for different down-type models; p and rj 
are the Wolfenstein parameters |^. 


From Table one can read that, in models of types s and b, it turns out that 
Br((i^)(f —)■ he) > Br(rf^)(f —)■ hu) while in models s we have the more exotic relation 
Br(rf^)(t — he) < Br(rf^)(f —)■ hu). The correlated allowed ranges for the branching 
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Figure 1: Correlations in flavour changing decays t ^ hu and t ^ he in models of type 
down quark - charged lepton, i.e. {di,ij). Line widths have no intrinsic meaning, they 
are only intended to help in distinguishing among the different models represented. 


ratios of these rare decay modes are shown in figure The correlations - the lines 
- are associated with each particular model; for example, the purple line labeled by 
the letters b and r is the correlation among the branching ratios Br(ft ,-)(t —)■ hu) and 
BY(^h,T)(t he), where the subscript (6, r) identifies completely the model - both the 
quark and the lepton type In that particular model, Br(fe {t —)■ he) can reach values 
up to a few times 10“^, that is up to the actual experimental upper bounds. In models 
(6, e) and (6, p) - in yellow and green respectively - the correlations among Br(t —)■ hu) 
and Br(t —)■ he) overlap along the same line with model {b, r), because all three models 
share the same V factors. The only difference between [b, r) and {b, /i) models in figure 
d , is in the maximum allowed value of the factor eisaitp +whose origin is in the 
aifferent predictions for the processes involving flavour conserving leptonic processes 
like h —)■ rr - as considered in appendix -. Without taking into account flavour 
changing low energy constraints, it is clear that the models (fe, r), (d, e), (d,/i), (s, e) 
and (s,/i), are the most promising models to discover new physics either in t ^ he or 
t —)■ hu. These models have flavour changing couplings in the up sector, therefore the 
Higgs coupling to uc could generate mixing 64 . To avoid too large 


mixing induced by tree level Higgs exchange, one can naively obtain the upper bounds 
for ej 3 a{ti 3 + collected in table 


Model d.y 

d 

s 

b 


< 0.05 

< 0.05 

< 72.3 


Table 5: Naive constraint on epa{tp + tjs^) from tree level h exchange. 
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The potential effects of these constraints are represented in hgure [T] with exclusion 
horizontal dashed lines for the corresponding models. The potential constraint in 6-type 
models is irrelevanlj^ while this constraint could be more relevant in s and d models. 
Nevertheless, in these models we do not have just the 125 GeV Higgs boson, but in 
addition another scalar H and a pseudoscalar A. It is well known that the scalar and 
pseudoscalar contributions to mixing cancel exactly in the limit of degenerate 

Note that the contributions to the oblique parameters 


masses 


66 do also cancel 


in the limit of degenerate masses and no mixing between the standard Higgs and the 
other neutral scalars 67 . These considerations imply that one cannot translate into 


Appendix B.2 


are not so evident, in Appendix B.3 


direct constraints the naive requirements on h couplings, since they can be relaxed in 
the presence of the other Higgses H and A. Although it is not within the scope of this 
paper to perform a complete analysis including the effects of the additional scalars, we 
illustrate how these cancellations operate in the case of meson mixing constraints, in 
and for the constraints coming from /i —)■ ey, where the cancellations 

In the following, potential low energy flavour 
changing constraints are shown in the figures in the same fashion as in figure It is 
important to keep in mind that they are indicative of which models are under pressure 
or else safer from that point of view. 

The correlations in figure correspond to the models of type down quark - charged 
lepton, {di,ij), where FCNC are present in the up quark sector and in the neutrino 
sector. In these models, FCNC are proportional to neutrino masses and thus there 
are no flavour changing constraints coming from the leptonic sector. The constraints 
from the Higgs signals involve the diagonal couplings which do change and were taken 
into account as explained in appendix When down quark-neutrino type models are 
considered, + 1^^)\ is also constrained by /i —)■ ey, r —)■ ey, r —)■ /ry and other 

flavour violating processes. It can be shown that, in any z/j model, /r —)■ ey is the most 
constraining process as far as Cj 3 a{t /3 is concerned. We address in more detail the 

/i —)■ ey restrictions in section [O] and in appendix |B. 3 


5.2 


Correlations in Up-Charged lepton models: h 
h —y bd 


bs versus 


In order to have h ^ bq decays at a significant rate, we have to consider up-type 


models, Uk, where FCNC occur in the down sector. In this case, following Eqs. (28) 


(|29|), the correlations among the Higgs flavour changing decays to down quarks are 

BT^u,){h ^ bs + sb) lUJ 


given by 


BT(^Uk){h ^ bd + db) 


(34) 


The values of the ratio in Eq. (34) for the different up-type models are given in table 


the correlations are represented in hgure 

The correlations in hgure follow from the full data analysis in appendix in¬ 
cluding the necessary study of the h total decay width in these BGL models. We can 


®Only for br models, where the top quark decay constraints |c/ 3 Q,(t /3 +tp^)\ < 4.9, Eq. ( [^ , can be 
saturated, perturbative unitarity requirements may be relevant (see appendix B.l). 


13 



















Model Uk 

u 

c 

t 

IWI" 

IWP 

AO = 0.0535 

i) = 18.688 



Table 6: V factors entering Eq. (34). 



Figure 2: Correlations in flavour changing h decays to down quarks in models. 


see that in {t, r), (c, e) and (c, fi) models, the channel h ^ bs + bs can have a branch¬ 
ing ratio at the 10“^ level, while in (c, r), {t,e) and (t,/i) models this branching ratio 
can be near 10“^ values. It is also remarkable that in [u, e) and {u, fi) models, the 
branching ratio of h ^ bd + bd can also reach values not far from the 10“^ level. Top 
endpoints of the correlations follow from the general analysis, which allows maximal 
values \ci 3 a{tp + t^^)\ ~ 8. As in the case of t —)■ hg correlations, the dashed lines 
show the naive constraints one would get from h contributions alone to the low energy 
processes - K^, - D^, and B^ - 5°; explicit values are collected in 

table Examples of models where low energy contraints can be relevant are models 
(m, e) and (n, p). Once again, we must stress that the presence of other contributions 
in these BGL models can relax these low energy constraints. We include them here for 
the sake of completeness. 


Model Uk 

u 

c 

t 

(t/S ) 

< 0.43 

< 0.43 

< 0.60 


Table 7: Naive - K^, B^ - B^ or 5° - 5° constraint on from tree 

level h exchange. 


Note, however, that since the values in Table are near 1, one cannot go too far 
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above the dashed lines in hgure without taking into account perturbative unitarity 
(see appendix |B.1[). 


5.3 Correlations in neutrino models: 
t hq or h bq' 


h fjLT together with 


In neutrino models, we have flavour changing Higgs interactions in the leptonic sector 
giving rise to the interesting processes h —)■ The corresponding cou¬ 

plings are proportional to one of the lepton masses, therefore the transitions including 
a r lepton are at least more probable by a factor We will concentrate in 

these transitions, containing a r lepton, even if experimentally the fie channel is very 
interesting. These transitions are also proportional to c^q ,(^/3 + , like all Higgs 

induced flavour changing transitions in these models; therefore, in down-type models, 
we will have perfectly defined correlations between h —)• /ir and t —)■ hq] in up-type 
models, we will have correlations among h ^ fir and h — bq. 


At present, as already mentioned, evidence from the CMS collaboration 56 points 
towards the possible observation of the decay h —)■ fir, which would definitely be 
“beyond the SM physics”. These predictions could be checked by looking at the cor¬ 
relations with the channels t ^ hq for down-type BGL models, and with h ^ bq in 
up-type BGL models. 

In BGL models of (d.^, z/g.) type, the correlations t ^ hq versus h ^ fif + rji, 


following Eqs. (19)-(23), are controlled by 


Br(rf^)(t ^ hq) 


2.063 


Vt, 


Y 

^Br(,^)(h^/if + r/i) 

h 


(35) 


Notice that these correlations are fixed by GKM and PMNS matrix elements. Never¬ 
theless, there is also the ratio of the total width of the Higgs in BGL models versus the 
SM value. This ratio makes the correlation to depart from strict linearity depending 
on CjSa and tp. In figure we first show the plot where only the range of variation 
is displayed, that is the strict linear relation. In this plot, figure one can see the 
effect of the upper bound on Br(i,^)(h fif + rfl). In particular, in models (5, z^i), 
the maximum value that Br(fe j^^)(f —)■ he) may reach is a few times 10“^. This value 
is smaller than the maximum allowed value in (6, r) models, presented in figure As 
usual, we have included - with dashed lines - the naive constraints coming from the 
individual Higgs contribution to low energy flavour changing hadronic processes. 

In BGL models of {uk, Ha) type, we have similar expressions for the correlation 
among h ^ bq and h ^ fif + rJi decays. The corresponding plot is shown in figure 
1^ We can observe again the effects of the measurement in the h —)■ /xr channel. 
Nevertheless, as one can see, h ^ bs + sb branching ratios can still have values above 
the 10“^ level. 

Although Figures!^ and 1^ show h ^ fir decays, the values corresponding to h —>■ er 
decays follow from the PMNS factors in Table Notice that for h —)■ e/x decays, an 
additional suppression factor — 3.5 x 10“^ is involved. It is important to 

stress that perturbative unitarity will not impose any further constraint on Figures 
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Figure 3: Correlations in flavour changing t —>■ he vs. h —)■ /rr decays in models. 


and [^because Eq. (25) is at work. Several authors have noticed that fj. ^ e'y constrains 
very severely the coupling h ^ fie via the mass unsuppressed two-loop Barr-Zee di¬ 


agrams 


Since in BGL models all leptonic flavour changing Higgs couplings are 
fixed by U, masses and +1^^), it is clear that the p —>■ ey bound will translate 

into an important constraint on + that has to be incorporated to the global 

analysis. However, in these two-loop diagrams - as in the case of the different neutral 
meson mixings not only the Higgs h can be exchanged, but also the other scalar H 
and pseudoscalar A will enter with the possibility to produce destructive interference 
between the different contributions. 


As detailed in appendix B.3 


we have made a 

global analysis fully incorporating the different contributions to /i —)■ ey (in BGL mod¬ 
els) in two different cases: ( 1 ) with varying free values of the masses mn and ruA below 
1 TeV, and (2) by imposing mn — m a < 50 GeV and varying mn below 1 TeV: in both 
scenarios, oblique corrections can be maintained under control. To illustrate how these 
cancellations operate in /i — )■ ey, we represent the correlation among h ^ fif + Tfi and 
f —)■ he in models (s, z/ 3 ) and (s, i^i). In figure we show this correlation in the full 
analysis, first without including the /i —)■ ey constraint, figure 5(a), and then, in figure 


5(b)| when we introduce the p —)■ ey constraint as mentioned in scenario (1), that is 


with free values of mn and rriA below 1 TeV. As figure]^ shows, the region of variation 
of the correlation remains essentially the same, meaning that there are cancellations 
at work, implying that in this kind of 2HDM, one cannot forget about the additional 
Higgses in order to impose the low energy constraints. Considering instead scenario 
(2), i.e. taking mn — mA < 50 GeV and varying mn below 1 TeV, the corresponding 
plot is shown in figure 

It is then clear that if we include in the analysis the complete two-loop Bar-Zee 
contribution to p —?• ey, we can conclude from the observed changes and the actual 
level of precision, that effects are not yet relevant in the majority of BGL models. We 
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Figure 4: Correlations in flavour changing h —>■ 6s vs. h ^ /ir decays in {u^, Ua) 
models. 


have illustrated this result with down-type models, but the same happens in up-type 
models; therefore, h —)■ 6s correlations with 6, —)■ /ir remain essentially unchanged 
without the inclusion of the /r —>■ ey constraint. 
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Figure 5: Correlations in flavour changing t —)■ he vs. h —)■ /ir decays in (s, ui) and 
(s, z/3) models. The regions, darker to lighter, correspond to 1, 2 and 3 ct regions from 
the full analysis. As anticipated, they are not simple straight lines. Furthermore, blue 
and red regions here correpond to yellow and purple lines in hgure The Icr region 
reflects the effect of the CMS measurement, which is compatible with zero at the 2.da 
level. 



Figure 6: Correlations in flavour changing t —)■ he vs. h —)■ /ir decays in (s, ui) and 
(s, z/3) models with /i —>■ ey following scenario (2) (see text and caption of hgure [^. 
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6 Conclusions 


We analyse flavour changing scalar couplings in the framework of a class of two Higgs 
Doublet models where these couplings arise at tree level, but with their flavour structure 
entirely determined by the CKM and PMNS matrices. This very special structure of 
the scalar couplings is stable under the renormalization group, since it results from a 
discrete symmetry of the Lagrangian. The symmetry can be implemented in various 
ways, corresponding to a variety of BGL models. We pointed out that this class of 
models leads to New Physics with potential for being discovered at the LHC and/or 
at an ILC. We examine in detail rare top decays like t ^ hq {q = u, c) and leptonic 
flavour changing decays such as h ^ {i = /i,e), as well as hadronic flavour 

changing decays like h ^ hs and h —)■ bd. All these decays occur in the SM only at 
loop level and therefore are strongly suppressed. In BGL models, the flavour violating 
couplings occur at tree level, but some of the most dangerous couplings are suppressed 
by small CKM elements. We address the question whether there are regions in some 
of the BGL models where these couplings are such that may lead to the discovery of 
rare flavour violating processes at the LHC-13TeV. 

We also do a systematic study of the correlations among various observables which 
are an interesting distinctive feature of BGL models. In the search of these regions, we 
have taken into account the low energy restrictions on flavour violating processes as 
well as the stringent constraints on all SM processes associated to the Higgs production 
and subsequent decays in the various channels. 

As far as the low energy flavour constraints are concerned, we agree with other au¬ 
thors that these cannot be imposed by assuming the dominance of the lighter Higgs con¬ 
tribution. This was known, in particular, in BGL models for neutral Meson-Antimeson 
mixing: there are important cancellations among different virtual Higgs contributions. 
We have illustrated this point showing how these cancellations operate in the two- 
loop, Higgs mediated, p —)■ ey process, where these cancellations can appear in the 
amplitude, operating at the level of one or two orders of magnitude. 

Two Higgs doublet models are one of the simplest extensions of the the SM. In 
general, they have a large number of free parameters and lead to scalar FGNG which 
have severe restrictions from low energy flavour violating processes. BGL models have 
the notable feature of having a small number of free parameters and achieving a natural 
suppression of these couplings, while at the same time allowing for the exciting scenario 
of having some flavour violating top and Higgs decays to occur at discovery level at 
the LHG-13TeV. 
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A Higgs signals 


Besides the appearance of flavour changing couplings of the Higgs boson, as shown in 


equations Eqs. (12)-(13), flavour conserving couplings are modihed owing to the mixing 


in the scalar sector, and thus a detailed analysis of the constraints on a — (3 and tan /? 
that Higgs data impose is mandatory. The experimental information concerning the 
SM-like with mass 125 GeV discovered at the LHC is summarised in a set of signal 
strengths of the form 


X _ [(Tij>V ^)]i 
hi ~ 


Br(h ^ X) 


a{pp -)■ h)sM\i Br(/i -)■ X) 


SM 


(36) 


where i labels the different combinations of production mechanisms and X the decay 
channels. Concerning production, the most relevant modes are gluon-gluon fusion 
(ggF), vector boson fusion (VBF) and Higgstrahlung (WH & ZH); values used in the 
analysis are collected in Table Relevant branching ratios within the SM are in turn 


Production channel i 

ggF 

VBF 

Wh 

Zh 

tth 

bbh 


a{pp -)■ h)sM, 

. (Pb) 

19.27 

1.578 

0.7046 

0.4153 

0.1293 

0.2035 


Table 8: Production cross sections (uncertainties are not shown), for rrih = 125.0 GeV 
and y/s = 8 TeV. 


collected in Table [H 


Decay channel X 

bb 

WW* 

ZZ* 

TT 

77 

99 

[Br(h —)■ X)]sM 

0.578 

0.216 

0.0267 

0.0637 

0.0023 

0.0856 


Table 9: Decay branching ratios for mh = 125.0 GeV (the total width is T^ = 4.03 
MeV). 


We now list the different signal strengths obtained by the CMS and ATLAS col¬ 
laborations organising them by decay channel. 


• h — 77 , 




1.12 


+0.37 

-0.32 


hvBF 


1.58 


+0.77 

- 0.68 


hvH 


77 

hggF 


1 'iO +0.41 
-0.33 


hvBF 


0.78 


+0.80 

-0.64 


-0 16 
U-J-O -0.79 


CMS pra 


ATLAS 71 . 


(37) 

(38) 


h-t ZZ 72 


5 


ho/ljet 


0.858 


+0.321 

-0.258 




1.235 


+0.852 

-0.583 


(39) 


where the dominant production in the 0/1 jet signal is gluon-gluon fusion - 
ggF:VBF ~ 20:1 -, while in the 2 jets signal, although gluon-gluon fusion con¬ 
tributes the most, VBF and VH are not negligible, - ggF:VBF:VH ~ 4:2:1 -. 
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h^WW 73 


,ww 

A^O/ljet 


0.74 


+0.22 

- 0.20 


,WW 

Mvef+vh 


0.60 


+0.54 

-0.46 


(40) 


where, as in /i^^, the 0/1 jet signal is gluon-gluon fusion dominated. 


h —)• rr 74 


/xSlet = 0.34 zh 1.09 , fill, = 1-07 ± 0.46 , fil]et-YBF = 0.94 ± 0.41, (41) 


where the 0 and 1 jet signals are dominated by gluon-gluon fusion and for the 2 
jet-VBF tag, ggF and VBF production are similar. 


h^bb 75,76 


= 1-0 ± 0.5 , /i^^H = 0.65 / 


bb 


, +0.43 


0.40 ) 


(42) 


Concerning the dependence of the couplings involved in the different production and 
decay channels, hWW and hZZ are rescaled by a factor with respect to the SM 
for all the models. This affects VBF and VH production modes, h —)■ WW, ZZ decays 
and the IF-loop contribution to h —)■ 77 , that we address later. For the couplings of 
h to fermions, the picture is more involved, they are modihed in a model dependent 

the changes in hit, hbb and hfr with 
respect to the SM, where the hff interaction is simply ^hff. 


manner. We remind in Tables 10 


11 


and 


12 


Model type 

u, c 

t 

Down model di 

Coupling rrit/v X 


^ pOL 

Sfia + C^a[(l Wtdi 


^0 


Table 10: htt coupling. 


Model type 

d, s 

b 

Up model u. 

Coupling 171},/V X 

^/ 3 q C/Sat/S 

^I3a 

S/3a + C/3o[(l \Vuib 

to 

1 

^Uib 

2f-il 


Table 11: hbb coupling. 


Model type 

e,/i 

T 

Neutrino model Ui 

Coupling rUr/v x 



S0a + Cpa\{\ jV/j 


^0 


Table 12: hrr coupling 


Notice that 

• the change in hfr only affects the branching ratio in 

• the change in hbb would in principle only affect the branching ratio in /if^; how¬ 
ever, production through the otherwise negligible bb ^ h process could be tan /3 
or tan“^ /3 enhanced. 
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• the change in htt affects gluon-gluon production and the top loop in h —>■ 77 
decays. 

The h —)■ 77 decay deserves some attention. In the SM it is a loop induced pro¬ 
cess where virtual W and top diagrams interfere destructively. Besides the individual 
rescaling of both contributions, additional contributions mediated by the charged scalar 
could also contribute. Scenarios with sizable contributions to h —)■ 77 require 
a specific analysis that is beyond the scope of this work. A regime with heavy 
bosons can always be considered where this approximation is sound. 

With all the ingredients in place, namely (i) the experimental constraints and (ii) 
the model predictions (simply expressed in terms of the different rescalings of SM 
couplings), a standard analysis of the {tan/?, a — /?} parameter space can be built. As 
an illustration of the effect of imposing agreement with the set of constraints on flavour 
diagonal Higgs couplings, we show allowed regions in the logio(tan/?) vs. a — (3 plane 
for a few models in Figure Since the overall agreement of different signal strengths 
with the SM is good, the region around Cf^a = 0 is in all cases allowed. Depending then 
on the particular structure of the tan (3 dependences in the couplings, the a — /3 span of 
the allowed regions for large or small values of tan /? can be anticipated. In addition it 
should be noticed that, in some cases, the fluctuations departing from signal strengths 
equal to 1, can be in fact accommodated with a — /? 7 ^ 7 r/ 2 . 



(a) Model t U 2 




1 

a — /3 


(c) Model u e 



Figure 7: 68 %, 95% and 99% CL regions in vs a — /? for a small sample of models. 
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B Constraints 


In this appendix we address the constraints imposed on by (1) pertnrba- 

tivity requirements for the couplings in the scalar potential, ( 2 ) tree level contributions 
mediated by the three neutral scalars h, H, A, to the mixing amplitude M 12 in neutral 
meson systems and (3) two loop Barr-Zee contributions to the p —)■ ey decay branching 
ratio involving the flavour changing interactions of all three neutral scalars, impose on 
cpaiip + It is important to stress that while h —)■ pr, t —)■ he, hu or h ^ bs, bd, 
depend on Cpa{ti 3 +and no other unknown parameter related to the scalar sector, 
the constraints analysed in sections B.l , B.2 and B.3 of this appendix, do involve new 
parameters like the masses mn and 


B.l Perturbative unitarity 

Neutral scalar masses and mixings arise from the scalar potential of the model and 
are related to the dimensionless quartic couplings A*. Perturbativity requirements like 
Aj < 47 r, could have some impact on the allowed values for cpa{tp + tp^)- Following 
appendix D of (here A 5 = Ae = Ay = 0, and furthermore notice in addition that 
in all Aj are two times our corresponding A* in Eq. (|^), 

'^\spaCpa = V^[s2a{^2S^p — Aic|) + (A 3 + Xj^S2pC2a)\- (43) 

Since cpa{tp + we have 


Sf3C^ 


— Cpa{tp + tp^)spa — —y ["S2«(A2t;3 — Ait^^) + 2(A3 + A4 )c2q:]. (44) 

TTl / 


It is then clear that, for rriA V, having cpa{tp + t^^) r\-/ 0(1) does not challenge naive 
perturbativity requirements like A* < dvr. For much larger values of m^, however, the 
situation is more involved, and only a detailed analysis including all relevant parame¬ 
ters can gauge the precise extent of the constraints on cpaitp + tp^) as a function of 
additional parameters. This is beyond the scope of the present work. Further relations. 


similar to (43), but involving tuh and ruh instead of tua, lead to the same conclusion 
on the perturbativity requirements for the Aj’s versus the values of cpa{tp -h tp^)- 
It is important to stress, however, that the presence of other constraints overrules the 
potential role of imposing perturbative unitarity in the scalar potential: h —)■ /ir alone, 
in neutrino type models, already requires cpaitp + t^^) ^ 1 ; bounds on rare top decays 
t —)■ hc,hu in models of types b and s impose cpa{tp + t'^^) < 5; for the remaining 
models constraints from meson mixings (addressed in the following), can play a more 
relevant role. 


B.2 Mesons mixings 

The contribution to the meson mixing amplitudes M 12 in BGL models is, to an excellent 
approximation (up to terms of order ^ for up models, analogously up to 


24 








—, — terms for down models), 

’ mt ' ’ 


Mi2 oc ((^ + ij‘)^ 




(45) 


In the strict = 0 limit, for = mu {a. sitnation which does not clash with 
electroweak precision data), there is a complete cancellation at work. Departing from 
the cpa = 0 limit, in the three dimensional parameter space, there is 

now a two dimensional snbspace where the cancellation is still complete 


{cU O = M\ [c^/mi + (1 - c^)/Mr^) 


Close to that snbspace, considering the h contribntion alone, does not reflect the actnal 
meson mixing constraint on c^q,. With respect to that single h contribntion, a cancella¬ 
tion of one or even two orders of magnitnde is achievable. However, larger cancellations 
that would be necessary in some cases to produce interesting phenomenological conse¬ 
quences (for example, for (u, e) and (n, /i) models in £gure|^, may be less likely. Since 
a detailed analysis would involve additional parameters, niH and only would-be 
bounds from h alone are explicitely shown in the different figures as an indication of 
potential constraint in some models. Other models, like all type b models, are free from 
such constraints. 


B.3 /i ^ 67 

The radiative transition /i —)■ ey in 2HDM is typically suppressed at the one loop level 


by lepton masses. Two loop contributions of the Barr-Zee type 68 where instead of 


two, only one suppressed scalar-fermion-antifermion coupling is involved - thus reduc¬ 
ing the chirality flip suppression - can be dominant. The contribution from this class 


of diagrams in the context of 2HDM was addressed in 78 . In BGL models, this two 
loop prediction has the general structure 

Br(/i -)■ 07 ) 2 loop = 

^ (^)'‘(V + l^(Q)t - 5.77 X i 0 -»(y + | 7 l,g,p, ( 46 ) 

where the neutrino-type model dependence is encoded in the prefactor \UejU'^j\ (which 
takes approximate values 0.32, 0.30, 0.11 for j = 1, 2, 3), and is the amplitude 

g 

A(q) = CfiaSfSa [^(rrih) - + -Kt [clJ{zh) + sIJ{zh) - g{,ZA)\ , (47) 

with yx = M^/Mx and zx = = h,,H,A), and Kt the model dependent 

change in the coupling of scalars to the top quark in the loop in Table IB ^ The 
functions /, g and h follow the definitions in [78], and 


E(m) = 3 /(s)+69(j/)+? |9(j/) + ft(!/)| + hh-h£) _ _ 


y = 


m? 


ml 


z = 




(48) 
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Model type 

u, c 

t 

Down dj 

Factor Kt 

tfi 


{^-\VttV)td-\Vtt\%^ 


Table 13: htt coupling, Kt factor. 


With the experimental bound a simple estimate only considering the h-mediated 
contributions and T,{mh) = 2.016 (neglecting the term Cp^f{zh) ioi a — j3 ~ 7r/2), puts 
the stringent bound CjSait^ +< 0.031. Nevertheless, ignoring the effect of the H 
and A terms is not justified. In Figure]^ we plot g{m‘^/rn?‘) and S(m) as 

function of m, together with H{mh) as a reference]^ Attending to Eq.(47) and Figure 



Figure 8: g{rnl/m?) and S(m). 


sizable cancellations among the different contributions can easily be at work. A 
simple exercise can be helpful to illustrate such scenarios: hxing mu = 500 GeV and 
requiring c^aitp + > 1 (largely in excess of the bound that h alone would give), 

one can numerically scan the range of values of ruA for which the experimental bound 
is satisfied. For an initial scanned range [100; 1000] GeV, values of uia satisfying the 
requirements cover the full range. In other words, deriving bounds from h alone is 
too simplistic. To remedy that and incorporate nevertheless the effect of /i —?• ey, we 
follow two steps: (1) allowed vs a — /3 regions are first obtained from a scan of Mh 
and Ma values which fulhll the experimental bound on Br(/i —)■ ey) for each model, 
(2) the general analysis of vs a — /? in each model is then restricted to the previous 
regions. For the hrst step, two different regimes are explored, (a) mH and tua are 
independent, (b) tuh and tua are required to be similar, namely rriA € [0.8; 1.2] x mn- 
Although the outcome of this auxiliary analysis is a reduction of the allowed tp vs. 
a — regions, prospects for the different rare decays under consideration in u models 
are not signihcantly altered. 

^Notice that, owing to the hierarchy in the CKM matrix elements, as a very good approximation, 
in down models d, s, Kt = tp, while in models 6, Kt = 

^Although Yj{m) would apparently induce non-decoupling contributions for large masses, E(to) 
grows as Inm [^, it appears multiplied by a cpa factor and thus, taking into account the previous 
discussion on perturbativity, there is no “real” non-decoupling. 
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